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We study two-photon exchange (TPE) in the elastic electron-nucleon scattering at high in 
the framework of pQCD. The obtained TPE amplitude is of order a/ as with respect to Born 
approximation. Its shape and value are sensitive to the choice of nucleon wavefunction, thus study 
of TPE effects can provide important information about nucleon structure. With the wavefunctions 
based on QCD sum rules, TPE correction to the electron-proton cross section has negative sign, 
is almost linear in e and grows logarithmically with up to 7% at — 30 GeV'^. The results 
of existing "hadronic" calculations, taking into account just the nucleon intermediate state, can be 
" smoothly connected with pQCD result near ~ 3 GeV^. Above this point two methods disagree, 

. which implies that "hadronic" approach becomes inadequate at high . Other relevant observables, 

' such as electron/positron cross section ratio, are also discussed. 
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] The two-photon exchange (TPE) in electron-proton scattering is actively discussed over last several years. The 
■ impetus for this was initially given by the discovery of so-called Ge/Gm problem in the proton form factor measure- 
ments [l| . It was shown later that the discrepancy between Rosenbluth separation and polarization transfer methods 
i; can be at least partially eliminated after taking into account TPE effects [2]. Several experiments aimed at direct 
Oh' detection of TPE contribution to cross section, are proposed 0]. Non-zero single-spin asymmetry, which is induced by 
Q^, the imaginary part of TPE amplitude, was observed experimentally The role of TPE in determination of proton 
^ radius 0, parity- violating observables fo"] and in deep inelastic scattering 0] was also discussed. 
f~| Currently, the measurements of proton form factors at ^ 10 GeV^ are in progress [1] and other measurements 

in this region are proposed [1, • Clearly, these experiments call for the reliable estimate of TPE effects for high-Q^ 
I kinematics, which was one of the aims of the present work. At moderate Q^, the TPE amplitude was calculated 
^ I using nucleon and resonances as intermediate states (further called "hadronic" approach) [2ij lU, J3j ■ At high Q^, 
\^ ' however, a natural means for the description of any process involving hadrons, and in particular TPE, is perturbative 
\^ , quantum chromodynamics (pQCD). Surprisingly, we have found no direct pQCD calculation of TPE in the literature. 
' In Ref. TPE at high was investigated using the formalism of generalized parton distributions. The authors 
] doubt of pQCD applicability in the presently accessible kinematical region and thus use an alternative method. The 
values of TPE corrections obtained this way have opposite sign to the results of "hadronic" calculations. The authors 
also use an assumption that the most important diagrams are that in which both photons interact with the same 
OO . quark. It turns out that in pQCD approach the situation is reversed (see below. Sec. IIII B|) . 

' In the present paper we study TPE in the elastic electron-nucleon scattering at high in the framework of pQCD. 
We employ the method, which was used to calculate baryon form factors in Refs. [l^, li6|. In the adopted approach, 
a nucleon with momentum p is represented as three collinearly moving quarks with momenta x^p, where < xi < 1, 
Si=i Xi = 1. All quark masses and nucleon mass are neglected and thus = [xip)"^ = 0. The process amplitude has 



' the form 

M^{c^{Vi)\T{y,,x,mx,)) (1) 

where T is hard scattering amplitude at quark level (represented by appropriate Feynman diagrams), (/'(x;) and 4>{yi) 
are initial and final nucleon spin- flavor-coordinate wavefunctions (quark distribution amplitudes). The convolution 
with nucleon wavefunction implies a convolution of spinor indices and an integration over dxidx2dx^5{l — xi — X2 — x^) 
and similarly for yi. 

To obtain non-zero transition amplitude one must turn the momenta of all three quarks from initial to final direction. 
In one-photon exchange (Born) approximation one therefore needs at least two hard gluons to be exchanged between 
the quarks. It follows then that the amplitude scales as aa1/Q^ and the nucleon form factor — as a^/Q*. In the case 
of TPE the exchange of one gluon is sufficient and thus the leading-order pQCD contribution to the TPE amplitude 
should be ~ a^as/Q^. The ratio TPE/Born then will be not just a, as one may naively expect, but a/as, which is 
significantly larger and growing with . Thus the larger is, the more important TPE will be. 

The paper is organized as follows. In Sec. [Til the observables, affected by TPE, are discussed, Sec. IIIII describes all 
ingredients of calculation (hard scattering amplitude for one-photon and two-photon exchange and nucleon wavefunc- 
tions), numerical results are given in Sec. llVl and conclusions — in Sec. fVl 



2 



II. KINEMATICS AND OBSERVABLES 

The momenta of particles are denoted according to e(fc) + N{p) e(fc') + N{p'). The transferred momentum is 
q = p' — p' , — —q^ > and ly ~ {p + k)(p' + k'). The reduced cross section of elastic electron-proton scattering 
can be written as 

an = ^^Qm\' + e\g,f + J^^l^^al^ (2) 
where M is proton mass, e = 1 - 2[1 + ly'^ / Q'^ {AAP + which for > AP turns to 

and Gmj Qe Etnd are certain invariant amplitudes (see details in Ref . [T3| ) . In Born approximation Qe and Qm 
become usual electric and magnetic form factors and vanishes, that is 

Ge ^ Ge + 5Qe, Gm = Gm + 5Qm, Gz = ^^3, (4) 

where the prefix 6 indicates TPE contribution. The dominant part of the cross section at high comes from the 
generalized magnetic form factor Gm and can be written as 

Hence we should primarily study TPE contribution 5Gm- This quantity also defines the positron/electron elastic cross 
section ratio: 



R=^- = 

(7 



Gm - 5Gm 



Gm + 5G 



M 



1 



l-4Re^ (6) 
Gm 



The generalized electric form factor is suppressed in the cross section by /Q"^ <C 1. Therefore a consideration 
of TPE corrections to Ge makes little sense. 

The amplitude SGs in principle can be measured in polarization experiments. Namely, neglecting terms of order 
M'^ /Q^ , the polarization of final proton in recoil polarization method is purely longitudinal and equals 

2e^ „ «3 



«./,Vr^(^l-jf^Rcgij (7) 

where h is electron helicity. Thus a precise study of e-dependence of may give an access to SGs ■ 

III. AMPLITUDE CALCULATION 
A. One-photon exchange 

At first, we briefiy review the pQCD calculation of form factors jl^, [l^. The nucleon form factors, or the elastic 
electron-nucleon scattering in one-photon exchange approximation, is described in pQCD to leading order in as by 7 
diagrams (Fig. [T^). For example, the piece of amplitude, coming from the first of them, is 



M = -^{^j • 6 • (4/9) • (-g2) • X 



(8) 



/J,/ M laiyiP +y2P -X2phl3{xip+q)'^t,(»Ja<»Jl3:>, 

xiy2X?,V3{yi + 2/2)(l - xi)q'^ 



where u and u' are electron spinors, the overall minus sign is due to negative electron charge, is ith quark charge, 
6 = 3! is a symmetry factor due to possible permutations of quark lines, {—q^) comes from quark and nucleon spinors 
normalization, 



4/9= (iA^iA^^iA-^^iA*) 



(9) 
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(b) 



(c) 



(d) 



FIG. 1: pQCD diagrams for eA'^ eN: one-photon exchange (a), two-photon exchange, leading order (b), subleading or- 
der (c,d). 



is a color factor (here A" are Gell-Mann matrices and ( ) means averaging over totally antisymmetric color wavefunc- 
tion). In the last equation as well as in Eq. ([5]) we separate matrices, acting on different quarks, by a (8) sign. Thus 
in the expression for color factor ^A^^A'' acts on the first, ^A" on the second and ^A*" on the third quark. 

Adding up contributions from all diagrams and using the fact that spin-flavor-coordinatc wavcfunction is totally 
symmetric under quark interchange, we obtain 

M = -^u'-ff^u ■ U'-f^U ■ Guiq^) (10) 

where U and U' are initial and final nucleon spinors and 

62 ei d 

xivixsysil ~ xi)(l - ys) x2y2X3y3{l - xi){l - ys) x2y2X3y3{i - X3){1 - yi) 

where hi — ±1 are signs of q uark hclicities; the helicities of initial and final quarks should be equal. This is equivalent 
to well-known result 




B. Two-photon exchange 

For the case of TPE, there are only 4 distinct diagrams in the leading order (Fig. [Ud), in which photons are 
connected to different quarks. The diagrams in which both photons interact with the same quark (Fig. [TJ;), need one 
more gluon to turn all quarks' momenta and thus are subleading in ag. Moreover, the evaluation of such diagrams 
alone is inconsistent, since the contribution of the same order in as comes from one- gluon corrections to the leading 
diagrams (e.g. Fig. [TJi). 

One point needs to be clarified here. If we remove the electron line, the diagrams Fig. [T]D-d will represent Compton 
scattering of virtual photons on the nucleon (doubly virtual Compton scattering, VVCS). And vice versa, TPE can 
be viewed as a process in which the virtual photon, emitted by the electron, is scattered from the proton and then 
absorbed back by the electron. VVCS has an important qualitative difference from the well-studied real Compton 
scattering (RCS). Since the momentum r of the real photon satisfies — 0, it cannot alone turn quark's momentum: 
{xip — yip'Y ^ {) = . Therefore diagrams with the structure like (Fig. ^p) vanish for RCS, and the amplitude 
expansion begins with 0{a^g) terms (diagrams like Fig. [T}:,d). On contrary, VVCS photons may be highly virtual, 
diagrams Fig. [1)3 contribute, and leading terms in VVCS amplitude arc 0{as). Hence one cannot employ an analogy 
with RCS in the analysis of TPE (cf. Ref. [H). 

We write down the expression for the first diagram in Fig. [T]d, the rest are analogous. We have 

SM = (^)'^-6.(-2/3).(-,^)x (12) 



u'-ffj,{k + X2P - y2V')lyU -fajyip' + y3P' - X3p)jf_, (g) 7^ (g) 7a | 

{k + X2P - y2P'Y + i'^ X2y2xly3{xi+X3){yi+y3Yq'^ 
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where the color factor is —2/3 = ^ (g) 1). After some algebraic transformations and using wavefunction 

symmetry, we obtain the full TPE amplitude in the form 

SM = -^u'-f^^u ■ U'^.U ■ (^^Sg^ + g^JGu) (13) 



q 



where 



fxr< xr \ 256Tr^aas , , 6162(1-/11/13) 

{SGM,Sg3) = {(j){y,) — 14 



1 f V — V + 



-2v,2v] \(j){x,)) 



V{X2 - y2) - q^{X2 + 2/2 - 2X2^2) +iO \1~ X2 1-^2 

From this expression it is easy to see the crossing symmetry of TPE amplitudes SG]\f and SG^: both are j/-odd. 
The quantity SQm, associated with the cross section correction ([S]), equals [l3| 

SGm - SGm + sSGs (15) 

As implied by Eqs. ([5][7]), it is better to consider ratios 5Qm/Gm and SQ^/Gm than the amplitudes themselves. This 
way we also avoid the uncertainty related with the absolute normalization of nucleon wavefunctions, since it cancels 
in the ratio. We have 

'SGm_ ^ 3a {<j>{y^)\{TsGM,T5g3)\(l){xi)) ^^^^ 



Gm'GmJ as {<p{yi)\TG,^,\(t>ixi)) 

where Tq^^^, Tsgm Tsg.^ are the expressions, sandwiched between wavefunctions in Eqs. (jllll4p . 

The obtained TPE amplitudes are free from infra-red (IR) divergence. This becomes clear if we recall that the 
IR-divergent terms are proportional to Born amplitude. Thus SQs is IR- finite (it vanishes in Born approximation), 
and 

SG^^^ ^ aGM\nX^ (17) 

where A is infinitesimal photon mass. The magnetic form factor [Eq. pT|) ] is a quantity of order O(a^). On the other 
hand, leading-order contribution to TPE is 0{as) and therefore IR divergence should only appear as a subleading 
effect, in the next order in as. 

Another interesting point pertains to photons' virtualities. In all diagrams they are both of order Q^, e.g. in the 
first diagram qf = —X2y2Q^ and g| — ~{xi + x-i){yi y3)Q^- We may conclude that the leading contribution to the 
amplitude at high comes from the region where both photons are hard, ql ^ q2 ^ Q^- 



C. Wavefunctions 



Before turning to numerical calculations, we must specify a model for wavefunctions. The requirement for total 
spin and isospin to be 1/2 together with Pauli principle fix the following form of the quark distribution amplitude 
(for the proton of positive helicity): 

In 

\(l}{xt)) = —^(l)i{xi,X2, X3) (lutuxdt) - \u^diu^)) + perm. (18) 

where "perm." means sum over all quarks permutations and /n is overall normalization constant not needed for our 
calculation. The neutron wavefunction is obtained by interchange d <-> u. As we can see, the distribution amplitude 
is completely determined by the function 0i . 

In general, the distribution amplitude and thus (pi depend logarithmically on Q^. Namely, we have 

MX^,Q^) ^ XiX2X3Y,['^'^(Q^)V'' BkPkixi,X3) (19) 
k 

where Pk are Appell polynomials (Pi = 1, P2 — xi ~ X3, etc.) and 7^ are corresponding anomalous dimensions 
p^ . Thus in the formal limit — > 00 the term with lowest 7^, which is Pi, dominates and (pi — > (pas = X1X2X3. 
This asymptotical wavefunction, however, leads to predictions inconsistent with experiment. In particular, it yields 
zero proton and positive neutron magnetic form factors. Thus at present the distribution amplitude should be 
considerably different from its asymptotic form [l^. Since the evolution with is very slow (7^: <C 1), the same 
wavefunction can be employed for all currently accessible with a reasonable accuracy. 

Various forms of distribution amplitude were proposed in the literature. We done the calculations with the following 
amphtudes: CZ 16], KS [17], COZ jS, GS [ig] and Het The CZ and KS amphtudes give practically the same 
results as the COZ amplitude, thus they are not considered further. 
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FIG. 2; TPE amplitude SGm/Gm vs. e at = 10 GeV^ 



FIG. 3: TPE amplitude SQm/Gm vs. e for neutron at 
= 5 GeV^ 




FIG. 4: Positron/electron cross section ratio for = 2,5, 10 and 20 GeV'^ (shown near the curves). Data are from Ref. (211 ] 
at 1.5 < < 5 GeV^ 



IV. NUMERICAL RESULTS 



There are two independent kinematical variables in any clastic process. For eiV scattering and e [Eq. ([3])] are 
generally used. The £-dependence of the obtained TPE amplitude 5Qm is shown in Fig.O It turns out to be universal 
for all (except for slow logarithmic evolution, which we neglect here). We see that the amplitude 5Qm^ calculated 
with COZ and Het wavefunctions, is very close to the linear function of e. Slight deviations from linearity are present 
near the cndpoints e — e — \ only. In contrast, GS wavefunction yields much larger and highly nonlinear TPE 
amplitude. In light of this it is worth noting that linear e-dependence of 5Qm is necessary and sufficient for Rosenbluth 
plots to remain linear even under the influence of TPE . Since careful analysis of experimental data do not reveal 
any nonlinearity in Rosenbluth plots [1^, we conclude that the experiment disfavors GS wavefunction. 

For the neutron target, both GS and Het wavefunctions yield nonlinear and anomalously huge TPE corrections, up 
to 25% (Fig. [3]). Taking into account the smallness of neutron electric form factor, these corrections would manifest as 
severe nonlinearities of Rosenbluth plots, that is, strong e-dependence of the elastic cross section. Though such cross 
section behaviour seems unlikely, the high-Q^ neutron form factor data are too poor to draw a final conclusion. Further 
experimental study of electron-neutron elastic scattering at high and different e can show definitely whether the 
nucleon is described by Het or by COZ wavefunction. For the present moment we take the COZ wavefunction as the 
most plausible. 

The amplitude SQ^, which determines the correction to longitudinal recoil polarization [Eq. ([7|], is small (< 1%) 
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FIG. 5: TPE amplitude SGm vs. at e = 0.5 (left) and £ = 0.1 (right). Dashed curves show hadronic calculations, with form 
factors parameterizations: dipole (red) and Ref. [24|] (black). 



for both proton and neutron, and unfortunately lies below the precision of today's experiments. 

The positron/electron cross section ratio is shown in Fig. [4l The calculation is done with COZ wavefunction at 
= 2, 5, 10 and 20 GeV^. The experimental data in the range 

1.5 < Q2 < 5 GeV^ from Ref. [2l| are also shown. 
Though the data points are well near the curves, the errors are very large. More precise data would be helpful, 
preferably in the low-£ region, where the predicted ratio is higher. 

The dependence of "normalized" TPE amplitudes at fixed e is completely determined by evolution of strong 
coupling constant as- We have used simple parameterization 

47r 

with A = 0.2 GeV. The resulting shape of TPE amplitude 5Qm for proton, calculated with COZ wavefunction, is 
plotted in Fig. [5] At « 30 GeV^, which is today the maximal ever investigated, the relative value of TPE 
amplitude reaches 3.5%, which corresponds to cross section correction of about 7%. Such a correction is however 
smaller than the errors of available data. On the other hand, TPE can be seen in recently proposed high-Q^ JLab 
experiment 0], where the estimated errors are at 1% level. 

The results of "hadronic" calculation [O, are also shown in Fig. [5] for comparison. Probably, the amplitude 
undergoes some gradual transition from this curve at lower to pQCD prediction at higher (recall that e- 
dependence in both cases is the same, approximately linear with positive slope). The figure suggests that a reasonable 
interpolation is possible between the "hadronic" result for below ~ 3 GeV^ and pQCD result above this value. 
But we also see a strong disagreement of these two curves at higher Q^. The most likely reason for such behaviour is 
that the "hadronic" approach, i.e. saturation of the intermediate hadronic states by the bare nucleon and the lowest 
resonances, is inadequate at high Q^. The multi-particle intermediate states yield a substantial part of the amplitude. 



V. CONCLUSIONS 



We have considered TPE for the elastic electron-nucleon scattering in the framework of pQCD. The calculations are 
done in the leading order with several model wavefunctions. For the proton target and wavefunctions based on QCD 
sum rules (CZ [Tf, KS [TtI and COZ Ts*]), the TPE amplitude 6Qm, which determines cross section correction, has 
linear e-dependence. Its value is of order a/as, grows logarithmically with and at = 30 GeV^ reaches 3.5% of 
Born amplitude. At lower a smooth connection is possible with previous "hadronic" calculations, in which TPE 
amplitudes were calculated taking into account just the nucleon intermediate state TT]. On the other hand, at high 
the results of these two methods are very different, which implies that "hadronic" approach becomes inadequate 
at Q2 > 3 GeV^ 

The size and e-dependence of TPE amplitudes are sensitive to the choice of nucleon wavefunction (quark distri- 
bution amplitude). At the same time, they are directly measurable: SQm/Gm via cross section or positron/electron 
cross section ratio and SQs/Gm — via longitudinal recoil polarization. Thus an accurate measurement of TPE ob- 
servables opens a new efficient way to study quark distribution amplitude in the nucleon. For example, the existing 
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experimental data already rule out GS wavefunction (Ref. fio']). Since TPE amplitudes have non-trivial e-dependence, 
they potentially provide much more information, than just nucleon form factors. Thus TPE turns from the correction 
to form factor measurements into an independent tool for studying nucleon structure. 
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